Memory effects on the statistics of fragmentation 
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We investigate through extensive molecular dynamics simulations the fragmentation process of 
two-dimensional Lennard- Jones systems. After thermalization, the fragmentation is initiated by a 
sudden increment to the radial component of the particles' velocities. We study the effect of tem- 
perature of the thermalized system as well as the influence of the impact energy of the "explosion" 
event on the statistics of mass fragments. Our results indicate that the cumulative distribution of 
fragments follows the scaling ansatz F(m) oc m~ a exp [— (m/mo) 1 ], where m is the mass, mo and 7 
are cutoff parameters, and a is a scaling exponent that is dependent on the temperature. More pre- 
cisely, we show clear evidence that there is a characteristic scaling exponent a for each macroscopic 
phase of the thermalized system, i.e., that the non-universal behavior of the fragmentation process 
is dictated by the state of the system before it breaks down. 

PACS numbers: 46.50.+a, 62.20.Mk, 64.60.-i 



I. INTRODUCTION 

The process of breaking solids into smaller pieces has 
been the subject of deep thoughts since the time of the 
Greeks, who tried to understand the building blocks of 
matter. Not going so far away in time or neither in the 
area of particle physics, the fragmentation process is still 
an important problem to study since it is a main issue 
in current problems in our day-to-day life. For instance, 
to understand why or how a material breaks is relevant 
in the development of new technological devices or in 
geological problems Q, Q • Because it is such a significant 
issue, a large number of experiments in fragmentation 
have been performed in order to collect data of fractures 
in many types of materials and objects forms 0,0, IS El 
0, 0- The number of theoretical articles on this topic 
is no smaller. The main focus of recent studies in this 
field is based on molecular dynamics (MD) simulations, 
where the results show an ubiquitous scaling behavior in 
the distribution of the mass fragments, F(m) ~ m~ a , 
with the exponent a depending on the dimensionality 
and initial parameters of the system 0, H, @ • 

The aforementioned experimental and theoretical stud- 
ies have shown that the mass distribution belongs to 
the same universality class for large enough input en- 
ergies when the MD system breaks into smaller pieces 
PTE 01 IE 113 • However using a molecular dynamics 
approach Ching et al. 22] fragmented an object repre- 
sented as a set of particles interacting via 

the Lennard- Jones potential with the fracture process 
being induced by 

random initial velocities assigned to the particles. The 
resulting steady-state form found for the cumulative mass 
distribution displays a typical power-law region, with 
a non-universal exponent that increases with the total 
initial energy given to the system. The same behav- 
ior has been observed in experimental fragmentation of 
long glass rods [23| and duly interpreted as an indication 
that the fragmentation process is not a self-organizing 
phenomenon, contrary to the assumption of Oddershede 



et al. |j] s In contrast to the self-organized criticality 
paradigm where the power law behavior should appear 
without a control parameter, there is a interesting claim 
that, in impact fragmentation, criticality could be tuned 
at a nonzero impact energy . In this way, the fragment- 
size distribution should satisfy a scaling form similar to 
that of the cluster-size distribution of percolation clusters 
[To| , but belonging to another universality class [Til Il2]| . 
From the results of such numerical models 

it has been suggested that there exists a critical im- 
parted energy, below which the object to be fragmented 
is only damaged, and above which it breaks down into 
numerous smaller pieces. The transition between the 
damaged to the fragmented states behaves as a critical 
point, with the fragment size distribution displaying a 
scaling form similar to that described in percolation the- 
ory [l(| ■ The same dependence on the impact energy 
for the fragmented state has been found very recently in 
another numerical model for the fragmentation of a circu- 
lar disk by projectiles |14| as well as in the experimental 
fragmentation of shells Q • 

In the present work our aim is to investigate through 
molecular dynamics simulations the effect of different ini- 
tial conditions (e.g., temperature and impact energy) on 
the mass distribution of fragments generated after an 
"explosion" takes place. One of our goals is to show that 
the scaling behavior observed in the statistics of mass 
fragments is non-universal and that this non-universality 
has a direct correspondence with the state of the system 
prior to fragmentation process. In Section II we describe 
the details of the model and simulations. The results 
are shown in Section III, while the conclusions and some 
perspectives are presented in Section IV. 



II. MODEL 

The fragmentation model used here is based on the one 
described in Refs. The initial state of the object 

to be fragmented is a thermalized configuration gener- 
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ated through a standard molecular dynamics simulation 
in the microcanonical ensemble. The particles interact 
with each other through a 6-12 Lennard- Jones pair po- 
tential and the system is brought to the desired equilib- 
rium temperature by integration of Newton's equations 
of motion 0|- A neighbor-list method is applied and 
periodic boundary conditions are used in all directions. 
This allows us to simulate up to 10 5 particles for a single 
realization of the fragmentation system. The results are 
then taken from an average of fifty realizations (the direc- 
tion of the initial velocities for the particles are different 
for each sample) for a given set of initial conditions, as 
defined by the value of the temperature, particle density, 
and energy given to break the MD system apart. This 
"explosion energy" is specified through the parameter R, 
defined as the ratio between the initial kinetic energy and 
the initial potential energy of the particle, immediately 
after the velocities are set according to the equation be- 
low 

v i (0)=vf + C7r i (0), (1) 

where vf are the initial velocities and Ti(0) are the ini- 
tial positions of the particles, obtained in the thermal- 
ization stage. The second term in the above equation is 
responsible for an expansion process that is preceded by 
an explosive event. The proportionality constant C has 
units of inverse of time and gives a measure of the initial 
energy imparted to the object. From time zero onward, 
no energy is added to the system and the particles' po- 
sitions and velocities are now calculated considering free 
boundary conditions. As a result, the system expands 
and the particles are distributed among clusters (frag- 
ments) of different masses. Each particle is considered 
as a monomeric cluster with unitary mass. There are 
several definitions for a particle cluster ^^|. Here, two 
particles will belong to the same cluster if they are sep- 
arated by a distance smaller than an arbitrary cutoff, 
r c = 3<7. The fragments are then classified according to 
their mass m and counted to enable the calculation of 
the distributions n(m) and F(m), both normalized here 
by the total number of fragments. 

III. RESULTS 

Due to the fluctuations in n(m), it is usually more 
convenient to work with the cumulative form of the mass 
distribution defined as |4| 

/>oo 

F(m) = / n(m!)dm! . (2) 
J m 

In Fig.^we show the behavior of F(m) for fixed values of 
the system density p = 0.61 and the parameter R — 0.43, 
and two different values of the temperature T. As can 
be seen, both distributions display a region of power-law 
behavior at intermediate values of m followed by a typical 
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FIG. 1: The cumulative mass distribution of fragments F(m) 
for p — 0.61, R = 0.43, and two different values of temper- 
ature T — 0.1 (circles) and 0.7 (squares). For comparison 
the dashed line corresponds to f(m) ~ exp [-(m/mo) 7 ]. The 
inset shows the corresponding mass distribution n(m), with 
the solid line representing the best fit. 

cutoff due to finite size. From Rcf. 16], the following 
expression has been proposed to describe the behavior of 
F(m): 

F(m) -TO 1 - Q exph(m/mo) 7 ], (3) 

where a is a scaling exponent, and mo and 7 are cut- 
off parameters. As depicted in Fig. ^ the application 
of a standard non-linear estimation algorithm to both 
data sets shows that Eq. © fits well the scaling region 
for intermediate masses as well as the decaying cutoff for 
large fragment sizes, which is compatible with a stretched 
exponential behavior. This fitting procedure gives esti- 
mates for the scaling exponent that are substantially dif- 
ferent, namely, a = 1.02 ± 0.01 and a = 1.20 ± 0.01 for 
T = 0.1 and 0.7, respectively. There is however a discrep- 
ancy between the data and the curve fitted with Eq. © 
at T — 0.7 for the region of small fragments. This can 
be readily explained in terms of the large "evaporation" 
rates at high values of the temperature - an expected ef- 
fect that is responsible for the progressive detachment of 
small clusters from the hull of large and medium frag- 
ments after the explosion event. The inset of Fig. 
shows that the corresponding behavior of the distribu- 
tion n(m) = dF(m)/dm for both values of T is also con- 
sistent with scaling ansatz Eq. J3J, when the parameter 
used are the same as those obtained for fitting its integral 
form F(m). 

In Fig. 12 we show the profiles of the distribution F(m) 
for several values of the temperature in the range 0.1 < 
T < 0.7 and a fixed value of p = 0.61. When observed 
in detail, the diversity in shape of F(m) for intermediate 
and large fragment sizes indicate that the fragmentation 
process must be restricted to a discrete and much smaller 
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FIG. 2: Log-log plot of F(m) for p = 0.61, R = 0.43, and 
different values of temperature. 




FIG. 4: Data collapse of the distribution F(m) shown in FigEl 
The collapse has been obtained by rescaling the abscissas m of 
each distribution to its corresponding cutoff parameters mo, 
as well as rescaling the values F(m) to F(mo). 
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FIG. 3: The behavior of the parameter mo against tempera- 
ture. The values of density and the parameter R are the same 
as in the previous figure. 



number of different classes of behavior than the variation 
with an entire spectrum of thermalization temperatures 
could suggest. This fact is quantitatively verified when 
we observe that, after fitting Eq. © to each data set, the 
scaling parameter a can assume one among only three 
distinct numerical values for the distributions generated 
at nine different temperatures. In Fig. [3] we show the 
variation of the cutoff parameter tuq with temperature 
for p = 0.61. As the temperature increases from T = 0.1, 
mo remains approximately constant up to T w 0.375, 
where it suddenly drops to become again constant, at 
least up to the maximum value of the temperature we 
use in our simulations, T = 0.7. This sharp transition 
in mo indicates the existence of a "critical" temperature 
below which a large cluster (i.e., with a size of the order 
of the system size) can exist. 

In Fig. 01 we show the data collapse obtained by rescal- 
ing the abscissas m of each curve shown in Fig. [3] to its 



corresponding estimate of the cutoff parameter mo, as 
well as rescaling the values F(m) to F(mo). This results 
clearly reveal the presence of only three groups of dis- 
tributions, each one with a characteristic value for the 
scaling exponent a. Such a behavior can be explained 
with the help of the phase diagram shown in Fig. [SI where 
the points following the vertical dashed line represent the 
values of temperature used in our simulations. From this 
diagram, we readily deduce that the threefold statistics 
of mass fragments shown in Fig.0|is a direct consequence 
of the three distinct phases to which the thermalized ob- 
jects belonged before they have been broken apart. It is 
interesting to note that, although the collapses are rather 
convincing for intermediate and large fragment sizes, the 
apparent divergence characterizing the statistics of small 
fragments due to "evaporation" appears to be continu- 
ously changing with temperature within each of the three 
groups of collapsed data. 

The situation becomes entirely different when we ana- 
lyze the influence of the energy parameter R on the statis- 
tics of the fragmentation process. In Fig. we show the 
distributions F(m) computed for MD systems thermal- 
ized with temperature T — 0.1, particle density p = 0.61, 
and for different values of R = 0.1, 0.2, 0.3, 0.4, 0.5, 
0.6, 1.0, 2.0, 3.0, and 4.0. From the non-linear fitting of 
Eq. 13 to each data set we notice that, while the scaling 
exponent remains approximately constant at a w 1.02, 
mo changes significantly with R. Precisely, as shown in 
Fig. {7\ the decay of mo with R can be described in terms 
of a power-law 



mo = a(R — Rq) 



-13 



(4) 



where a — 640.0±0.1 is a prefactor and the exponent /? = 
0.67 ±0.02. The parameter Ro is an offset that is related 
to the competition between the thermal energy of motion 
and the energy that holds the system together, i.e., the 
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FIG. 5: The phase diagram for a two-dimensional system 
with particles interacting through the Lennard- Jones poten- 
tial. The points following the vertical dashed line represent 
the values of temperature used in our simulations. 
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FIG. 7: Log-log plot of the crossover parameter mo against 
the energy difference (R — Ro) for MD systems thermalized 
with p = 0.61 and T = 0.1. 




FIG. 6: Log-log plot of the cumulative distribution F(m) for 
values of the energy parameter R ranging from R — 0.1 to 
4.0. The MD systems have been thermalized with T — 0.1 
and p = 0.61 



ratio between the kinetic energy and the potential energy 
just before the velocities are settled according to Eq. Q] 
and the boundary is lifted. Using Eq.^Jand its estimated 
parameters to rescale the data presented in Fig. we 
show in Fig. |H| that the distributions for all values of 
R can be nicely represented by a single data-collapsed 
curve. Of course, this should only be valid for systems 
subjected to the same thermalization process, i.e., if p 
and T are kept constant. 



IV. CONCLUSIONS 



ular dynamics simulations. Specifically, we have shown 
how the statistics of the fragmentation process depends 
on (i) the thermalization temperature of the system be- 
fore its breakdown, and (ii) the energy imparted to the 
system to induce fragmentation. In the first case, we 
verified that the cumulative mass distribution follows a 
power law for intermediate masses, with an exponent that 
depends on the region of temperature considered. More 
precisely, we showed that it is the phase of the thermal- 
ized object that is responsible for the difference in these 
scaling exponents. It means that the process studied here 
can be rather sensitive to the previous state of the sys- 
tem, although it introduces a significant disturbance from 
an energetic point of view. As a consequence, fragmen- 
tation carries memory. Finally, we turned our attention 
to the variation of the parameter R. Differently to the 
previous case, we obtained a unique scaling exponent for 
the cumulative mass distribution for different values of R. 
This result is in good agreement with previous studies in 
the literature indicating some sort of universal behavior 
present in fragmentation processes |bl| . 
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FIG. 8: Data collapse of the distribution F(m) for p — 
0.61, T — 0.1, and different values of the energy input 
R. For comparison the dashed line corresponds to f(m) ~ 
exp [-(m/m ) 7 ]. 



[1] D. L. Turcotte, J. of Geophys. Res. 91 B2, 1921 (1986). 
[2] F. Wittel, F. Kun, H. J. Herrmann, and B. H. Kroplin, 

Phys. Rev. Lett. 93, 035504 (2004). 
[3] T. Ishii and M. Matsushita, J. Phys. Soc. Jpn. 61 3474, 

(1992). 

[4] L. Oddershede, P. Dimon, and J. Bohr, Phys. Rev. Lett. 

71, 3107 (1993). 
[5] A. Meibom and I. Balslev, Phys. Rev. Lett. 76, 2492 

(1996). 

[6] M. Marsili, and Y.-C. Zhang, Phys. Rev. Lett. 77, 3577 
(1996). 

[7] T. Kadono, Phys. Rev. Lett. 78, 1444 (1997). 

[8] H. Inaoka, E. Toyosawa, and H. Takayasu, Phys. Rev. 

Lett. 78, 3455 (1997). 
[9] F. Kun and H. J. Herrmann, Phys. Rev. E 59, 2623 

(1999). 

[10] D. Stauffer and A. Aharony, Introduction to Percolation 
Theory (Taylor Francis, Philadelfia, 1994). 

[11] J. A. Astrom, B. L. Holian, and J. Timonen, Phys. Rev. 
Lett. 84, 3061 (2000). 

[12] H. Katsuragi, D. Sugino, and H. Honjo, Phys. Rev. E 68, 
046105 (2003). 

[13] U. Naftaly, M. Schwartz, A. Aharony, and D. Stauffer, J. 



Phys. A: Math. Gen. 24, L1175 (1991). 

[14] B. Behera, F. Kun, S. McNamara, and H. J. Herrmann, 
|arXiv:cond-mat /0404057| 

[15] A. Diehl, H. A. Carmona, L. E. Araripe, J. S. Andrade 
Jr., and G. A. Farias, Phys. Rev. E 62, 4742 (2000). 

[16] A. Diehl, J. S. Andrade Jr., and G. A. Farias, Phys. Rev. 
E 65, 048102 (2002); J. A. Astrom, R. P. Linna, and J. 
Timonen, Phys. Rev. E 65, 048101 (2002). 

[17] D. C. Rapaport, The Art of Molecular Dynamics Simu- 
lation (Cambridge Press, United Kingdom, 1997). 

[18] N. Sator, Phys. Rep. 376, 1 (2003). 

[19] X. Campi, H. Krivine, E. Plagnol, and N. Sator, Phys. 

Rev. C 67, 044610 (2003). 
[20] F. Kun, and H. J. Herrmann, Phys. Rev. E 59, 2623 

(1999). 

[21] X. Campi, H. Krivine, and N. Sator, Physica A 296, 24 
(2001). 

[22] E.S.C. Ching, Y.Y. Yiu, and K.F. Lo, Physica A 265, 
119 (1999). 

[23] E. S. C. Ching, S. Liu, and K.-Q. Xia, Physica A 287, 
83 (2000). 



